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ABSTRACT 

Here we present previously unpublished optical spectra of supernova (SN) 2001ig, a Type lib SN, 
from about a week after explosion until nearly one year later. The earliest spectrum consists of only 
a few broad absorption features, but soon more common Type II SN features including hydrogen 
P-Cygni profiles and helium absorption become apparent. At later times, as the H features fade and 
the He I absorption becomes more prominent, we observe the SN to transition from a Type II to a 
Type lb. Finally, observations after 250 days past explosion show a nebular-phase SN spectrum with 
one of the largest magnesium to oxygen intensity ratios ever seen. Additionally, we present models 
of the late-time spectra which indicate that the inner ejecta consist of ~1.15 M Q of material, most of 
which (by mass) is in the form of oxygen, with ^0.13 M Q of 56 Ni and essentially no hydrogen. 

Subject headings: supernovae: general — supernovae: individual (SN 2001ig, SN 1987K, SN 1993J, 
SN 2008D, SN 2002ap) — radiative transfer: supernovae 



1. INTRODUCTION 

It is thought that most high-mass stars (> 8 M Q ) 
end their lives as core- collapse supernovae (SNe) (e.g., 
IWooslev fe Jankal 12005, and references therein). These 
SNe are divided spectroscopically into subgroups based 
on the strength of H and He in their optical spectra. 
The progression from Type II to lb to Ic corresponds to 
SNe having both H and He, to just He, to ne ither H nor 
He p resent in their spectra (for a review, sec Filippenko 
[1993). The class of Type lib SNe (SNe lib) have spectra 
which undergo a transition from Type II to Type lb as 
their H features fade with time. It is likely that SNe lib, 
lb, and Ic undergo significant mass loss, due either to 
stellar winds and eruptions or to mass transfer to a bi- 
nary companion, before explosion. 

SN 2001ig was discovered i n the nearb y late- type 
spiral galaxy NGC 7424 by lEvans et all ([20011 ) on 
2001 Dec. 10.43 (UT dates are used throughout this pa- 
per). It is loca ted at a.i2oop = 22 h 57 m 30f74 and <5j2ooo = 
-41°02'26'/1 (|Rvder et all 120011 ). A n explosion d ate o f 
2001 Dec. 3 has been estimated by iRvder et ail (|2004f ) 
from modeling the radio light curve, and we use this as 
the reference date for our spectral observations through- 
out this paper. A spectrum was obtained within two 
days of detection of the SN and showed only a few broad 
absorption features. However, the spect rum from Las 
Campanas Observatory on 2001 Dec. 13 (|Phillips et al.1 
12001ft revealed si milarities to the Type lib SNe 198 7K 
(Filip penkdll988l) and 1993J (|Filippenko et al.lll993ft . 

Over the following month, spectra from the Euro- 
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pean Southern O bservatory (|Clocchiatti fe Prietoll200ll 
IClocchiattil 12002ft showed SN 200 lig to evolve in a 
manner similar to both SNe lib 1987K and 1993J, 
as the hydrogen lines became weaker and showed the 
complex shape charact eristic of the transitional stage 
(jFilippenko et al.lll994f ). By Oct. 2002, the predicted 
transition to a Type lb SN in the neb ular phase was 
complete (Filippen ko fe Chornock! [2002) . At the time 
it was also noted that the Mg I] A4571 feature was the 
strongest (relative to [O I] AA6300, 6364) ever seen in a 
SN. Our data complement the optical spectropolarimetry 
an d total-flux spectra th at were presented and analyzed 
bv lMaund et~ai1 (|2007bl) . 

In addition to optical data, X-ray images were ob- 
tained with the Advanced CCD Imaging Spectrometer-S 
(ACIS-S) instrument on the Chandra X-ray Observatory 
on 2002 May 2 2, showing a 0.2-10.0 k eV luminosity of 
-10 38 ergs s" 1 (jSchlegel fc Rvdenl27)02ft . SN 2001ig was 
also detected at radio frequencies; f or further informatio n 
on the radio data and analysis, see lRvder et al.l (|2004ft . 

Here we present and analyze optical spectroscopic data 
as well as two spectral models for SN 2001ig at late times. 
In 32] we describe our observations and data reduction, 
and in $3] we discuss our analysis of the spectra. Our 
model spectra and their implications are presented in iJH 
We summarize our conclusions in ij5] 

2. OBSERVATIONS AND DATA REDUCTION 

Beginning one week after explosion, optic al spectra 
were obtained using the LDSS-2 spectrograph (jMulchae vi 
2001) mounted on the Magellan Baade 6.5-m tele- 
scope, the red arm of th e dual imager/spectrograph 
EMMI (Dckke r et al.l Il986) mounted on the European 
Southern Observatory (ESO) 3.6-m New Technology 
Telesco pe, the L ow-Resolution Imaging Spectrometer 
(LRIS; lOke et al.l 11995ft mounted on the Keck I 10-m 
telescope, and the d ual imager/spectrograph EFOSC2 
([Buzzoni et al.lll984ft mounted on the ESO 3.6-m tele- 
scope. Our last spectral observation occurred on 
2002 Nov. 8, approximately 340 days after explosion. All 
observations used a 0.7"-1.0" wide slit and were aligned 
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along the parallactic a ngle to reduce differential light 
losses (|Filippenkol fl982). Table Q] summarizes the spec- 
tral data of SN 2001ig presented in this paper. 

All spectra were red uced using standard techniques 
(e.g.- lFolev eTafl 12003). Standard CCD processing and 
spectrum extraction for the LRIS data were completed 
with IRAF 8 , and the da t a wer e extracted with the opti- 
mal algorithm of iHornel (|1986f ). We obtained the wave- 
length scale from low-order polynomial fits to calibration- 
lamp spectra. Small wavelength shifts were then applied 
to the LRIS data after cross-correlating a template sky 
to the night-sky lines that were extracted with the SN. 
Using our own IDL routines, we fit spectrophotomet- 
ry standard-star spectra to the LRIS data in order to 
flux calibrate our sp e ctra and to remove t elluric lines 
(|Wade fe Hornd[l988HMatheson et al.ll2000D . 

Each EMMI observation was reduced independently 
and the final spectra were combined with a pixel- 
rejection algorithm based on the mean. 

3. SPECTRAL ANALYSIS 
3.1. Early-Time Spectra 

We present our early-time spectra of SN 2001ig in Fig- 
ure [TJ all have had the rec ession velocity of NGC 7424 
removed (939 km s _1 , from Koribalski et al. 2004). The 
first spectrum of SN 2001ig, taken on 2001 Dec. 11 (8 
days past explosion) , shows only a few strong, broad ab- 
sorption features. We conjecture that the feature near 
4000 A is H7 and the feature near 4300 A is a blend of 
Fe II lines (near a rest wavelength of 4570 A) blueshiftcd 
by ~ 22,500 km s _1 and ~ 18,200 km s _1 , respectively. 
Unfortunately, this identification is somewhat dubious 
since the H7 feature appears to be unusually strong. 

Instead, we note that this pair of features resembles 
the "W" feature se en in very early-tim e spectra of the 
Type II SN 2005ap (lOuimby et al.|[2007h and the Type lb 
SN 2008D (|Malesam et all 120091 : iModiaz et all 120091 a 
spectrum of which is reproduced in Figure [TJ. The fea- 
tures in both objects were attributed to a combination 
of C III, N III, and O III at blueshifted velocities of 
-21,000 km s" 1 and -30,000 km s" 1 , respectively. If 
we assume that these are the correct line identifications, 
then they are blueshifted by —26,000 km s _1 in the day 
8 spectrum of SN 2001ig. However, there is some uncer- 
tainty in the bluer of these two features since its left wing 
is distorted somewhat by the noise at the bluest end of 
our sp ectrum. It should also be noted that lMazzali et al] 
(2008) found no evidence for the "W" feature in their 
earliest spectrum of SN 2008D. 

The third feature (near 4600 A) can be decomposed 
into two Gaussian profiles centered at — 4498 A and 
-4596 A, which could be a blend of H/3 and He I AA4921, 
5015 blueshifted by -23,200 km s" 1 and -23,300 km s _1 , 
respectively. At this epoch SN 2001ig defies sim ple clas- 
sification (as pointed out by Phillips et al. 200l|). 

There is also weak, but broad, absorption at —5300 A 
and 5420 A, the latter of which could be He I A5876 
blueshifted by — 24,100 km s -1 . In addition, the ab- 

8 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science 
Foundation. 



sorption near 6060 A is most likely Ha blueshifted by 
-24,000 km s -1 . 

We measure a blackbody-fit temperature of —14,000 K 
for the continuum of this observation. In comparison, the 
derived early-time blackbod y temperatures of SNe 1 993J 
and 2008 D were -12,000 K (iClocchiatti et al.l ll995h and 
13,000 K (jModiaz et alj|2009 ), respectively, which match 
nicely with SN 2001ig. 

In our highest-resolution spectrum, taken on 
2002 Jan. 18 (46 days after explosion; see H3.2I for 
more on this observation), we measure an upper limit 
to the equivalent width (EW) of the Na I D line in 
NGC 7424 of —0.1 A. This correspon ds to a host-galaxy 
reddending of E (B - V) < 0.02 mag (jMunari fe Zwitterl 
119971) . which is comparable t o the Galactic r edden ing of 
E (B — V) — 0.011 mag from lSchlegel et al.1 (|l99l . The 
low host reddening for SN 2001ig is reasonable since the 
SN occurred on the outskirts of the galaxy Since these 
numbers are both quite small, we ignore re ddening when 
plotting our spectra. We also note that iMaund et al.l 
(2007b) required a total reddening somewhat in excess 
of the Galactic value plus the value we calculate for the 
host to explain the polarization they observe. 

In all other early-time spectra, we detect He I A5876 
and P-Cygni profiles of Ha (with the absorption compo- 
nent greatly dominating over the emission component). 
The observation from 2001 Dec. 13, for example, shows 
the He I feature blueshifte d by —16,000 km s -1 . In addi- 
tion, iPhillips et al.l (|2001f ) mention that there are broad 
absorption features at 4080 A and 4580 A in the observa- 
tion from 2001 Dec. 13, and it is possible that they are H7 
and H/3, respectively. However, the observation lacked 
acceptable blue-side flatfields and proper wavelength- 
calibration files, so data at wavelengths <5000 A cannot 
be completely trusted. We also identify the very weak, 
broad absorption near 7300 A as O I A7774. 

By 13 days after the explosion many more spectral fea- 
tures have appeared. Still present are the P-Cygni pro - 
file of Ha and the He I A5876 line. IMaund et all (|2007bD 
point out that Ha at this time is likely to be blended 
with He I A6678. No longer visible in this spectrum is the 
aforementioned "W" feature that probably came from a 
combination of C III, N III, and O III. This is again sim- 
ilar to SN 2008D, whose "W" disappeared in the span 
of -1 day (jModiaz et al.ll2009h . However, there are still 
features at about 4130 A and 4287 A, and we associate 
these with possible absorption due to H7 and Fe II, re- 
spectively. 

The feature near 4600 A is most likely H/3 while the 
noisy dips just below 4800 A are probably due to a com- 
bination of He I AA4921, 5015 and Fe II AA4924, 5018. 
Also seen is weak absorption likely due to Fe II A5169 
near 5030 A and quite broad, weak absorption from the 
Ca II near-infrared (IR) triplet around 8150 A. 

On 2001 Dec. 23, we acquired the fifth spectrum from 
the top of Figure [TJ Again, we detect the Ha P-Cygni 
line, although now it has a far more complex profile with 
a larger full width at half-maximum intensity (FWHM) 
than previously seen. This implies that the Ha feature 
is now likely becoming blended with Si II A6355. The 
feature near 5600 A in this spectrum clearly shows two 
local minima which are probably due to He I A5876 and 
Na I D. Also in this observation we detect distinct H/3, 
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TABLE 1 

Journal of Observations 



UT Date 


Day a 


Telescope 13 


Range 


Res. 


Airmass d 


Exp. 








(A) 


(A) c 




(s) 


2001 Dec. 11 


8 


LDSS2 


3800-9000 


13 


1.07 


2 x 90 


2001 Dec. 13 


10 


LDSS2 


5000-10000 


13 


1.18 


2 x 60 


2001 Doc. 16 


13 


EMMI 


3980-9400 


9 


1.64 


3 x 60 


2001 Dec. 23 


20 


LRIS 


3400-10300 


7 


2.49 


60 


2002 Jan. 13 


11 


EFOSC2 


3360-10300 


9 


1.87 


3 x 90 


2002 Jan. 17 


15 


LRIS 


3520-10000 


7 


4.09 


360 


2002 Jan. 18 


46 


LRIS 


3500-5380 


3 


4.99 


450 


2002 Jan. 18 


16 


LRIS 


5710-7000 


2 


4.99 


150 


2002 Jan. 29 


57 


LDSS2 


3600-10000 


13 


2.61 


90 


2002 Aug. 16 e 


256 


FORS1 


4000-8600 


12 


1.04 


4 x 1200 


2002 Oct. 8 


309 


LRIS 


3150-9420 


6 


2.18 


2 x 600 


2002 Nov. 8 


340 


LRIS 


3100-9410 


6 


2.05 


900 



a Days since the explosion date of 2001 Dec. 3 (estimated by Rvdcr ct afl l2004l . 
from t he radio light cu rve). LDSS2 = Magellan Baade 6.5-m/LDSS-2 spectro- 
graph (Mulchacy 2001); EMMI = European Southern Observatory (ESO) 3.6- 
m New Technology Telescope/red arm of the dual imager /spectrograph EMMI 
dDekker et a l. 1986); LRIS = Keck I 10-m/Low- Resolution Imaging Spectrome- 
ter IIQke et a l. 1995); EFOSC2 = ESO 3.6-m/dual imager /spectrograph EFOSC2 
l IBuzzoni et al.ll984l) ; FORS1 = ESO 8.2-m Very Large Telescope/FORSl spectro- 
graph ( Appcnzcllcr ct al. 1998).° Approximate spectral res olution, A. d A i rmass a t 
midpoint of exposures. 6 Spectrum previously published by Maund ct al. (2007b). 
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Fig . 1. — Spectra of SN 2001ig (and three comparison SNe) with some line identifications. The top spectrum is SN 2008D (IModiaz et al.l 
2009) and has been dereddened by E (B - V) Host = 0.6 mag. The next four spectra, of SN 2001ig, from top to bottom were obtained 
on 2001 D ec. 11, 13, 16, and 23 ; ages shown ar e relative to the adopted explosion date of 2001 Dec. 3. The bottom tw o spectra are 
SN 1993J dFilippenko et alj|1994f ) and SN 1987K jFiIippenkdll988l) . Recession velocities of 2100 km s" 1 HModiaz et al.ll2009T) . 939 km s _1 
(Koribalski ct al. 2004), —34 km s _1 , and 805 km s — 1 (both from dc Vaucoulcurs et al. 1991) have been removed from the spectra of 
SNe 2008D, 2001ig, 1993J, and 1987K, respectively. Telluric features are indicated by the © symbol. 
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He I A4921, and He I A5015 absorption, in addition to 
absorption from the Ca II near- IR triplet around 8200 A. 
Finally, extremely strong absorption due to Ca II H&K 
can now be seen near 3750 A. 

It has been pointed out that in many ways the 

first few spectra of SN 2001ig resemble tho se of 

SNe 1987K and 1993J at ea rly times (|Phillips et al.ll2001l ; 
IClocchiatti fc Prietoll200il ). Examples of both at similar 
ages are plotted for comparison in Figure [1] and there do 
appear to be quite a few similarities among these objects. 
However, it should be noted that Ha (both the absorp- 
tion and emission components) is weaker in SN 2001ig. 
In addition, the overall continuum shape of SN 2001ig 
does not change much during this two-week period. 

Figure [T] also indicates that SN 200 lig has a bluer 
continuum than either SN 1987K or SN 1993J; how- 
ever, this could be due to reddening. Both SNe 1987K 
and 1993J are plotted with no reddening correction, 
which is approximate ly consistent with observations 
(|Richmond et al.lll994h . but there is also evidence that 
the color excess of SN 1993J is as high as 0.25 -0.32 mag 
(|Richmond et al.lll994HClocchiatti et al.lll995[) . 

3.2. Transition Spectra 

Our spectra from ~5.5 weeks past explosion to ^8 
weeks past explosion are presented in Figure [2j The most 
obvious change in the spectra of SN 2001ig from FigureQ] 
to Figure [2] is that the flux in the blue part of the con- 
tinuum has decreased significantly; SN 2001ig now looks 
very much like SN 1993J in its overall continuum shape. 
In addition, SNe 2001ig and 1993J look closely similar 
in the range ~ 4000-6000 A, much of which is a forest 
of Fe II blends (with som e hydrogen Balme r and He I 
lines as well). As noted bv IClocchiat ti (2002), the main 
difference at this time between SNe 2001ig and 1993J is 
that SN 2001ig has stronger Ca II H&K lines, weaker H I 
and He I lines, and f aster expans i on ve locities. 

Also mentioned bv lClocchiattll (|2002D . by 2002 Jan. 13 
(the top spectrum in Figure [J), SN 2001ig has begun its 
transition from Type II to Type lb. This transformation 
is marked by the decrease in strength of Ha and H/3 
while the He I spectral features increase in strength. The 
most notable cases of this are the dramatic appearances 
of the He I A6678 absorption feature, which appears as 
a notch on the red wing of the Ha absorption profile, 
and strong He I A7065 absorption, which unfortunately 
coincides with a telluric absorption band in the day 41 
and 57 spectra. This transformation continues during 
the subsequent two weeks of observations (i.e., the rest 
of the SN 200 lig spectra in Figure [2]). 

The minimum of the Ha absorption feature is be- 
coming more complex and appears to have developed a 
doubled-troughed profile. This is to be expected since at 
this time Ha is being heavily ble nded with Si I I A6355 . 
In addition, as was pointed out by lMaund et al.l (|2007bl ) . 
the He I A5876 feature has become more prominent and 
the red wing of its absorption now shows a notch due to 
Na I D absorption. 

Also, by this time the Ca II H&K feature has devel- 
oped a P-Cygni profile with quite a strong emission com- 
ponent. In addition, the Ca II near-IR triplet absorption 
has strengthened dramatically. The previously identified 
absorption from Fe II near 4300 A and 5000 A becomes 
stronger, broader, and more complex with time. This 



is likely du e to line blanketi ng which was also seen in 
SN 2002ap (|Folev et al.H2003T ). 

The broad absorption near 9000 A could be O I A9266 
at about —9000 km s -1 . In addition, strong, relatively 
broad absorption from O I A7774 is again seen in these 
observations, but at roughly —8000 km s _1 . However, 
this feature coincides with a telluric absorption band in 
the day 41 and 57 spectra, making the velocity determi- 
nation difficult and uncertain. 

3.3. Late-Time Spectra 

A late-time spectrum from lMaund et al . (2007b) (256 
days past explosion) and our late-time spectra (309 and 
340 days past explosion) are presented in Figure [3] and 
are compared with similar late-time spectra of SNe 1993J 
and 2002ap. SN 2001ig has reached the nebular phase 
by this time and has quite strong [O I] AA6300, 6364 
and Mg I] A4571 emission. These and other common 
nebular SN features are marked in Figure [3j the features 
of SNe 2001ig, 2002ap, and 1993J are all similar, but the 
relative strengths of some of the emission lines differ. 

One of the most notable differences among the ob- 
jects is that SN 1993J exhibits obvious Ha emission 
by day 298, appearing as a shoulder on t h e red wing 
of [O II AA6300, 636 4. iFilippenko et alj (|1994h and 
iMatheson et al.l l|2000f ) showed that the Ha feature grad- 
ually strengthened from this time onward in SN 1993J, 
whereas in spectra at similar epochs of SN 2001ig shown 
in Figure [3] there is no obvious Ha emission. For 
SN 1993J the hydrogen at late times was attributed to 
interaction with circumstellar gas, but in SN 2001ig we 
see no evidence for such an interaction. 

Figure [3] also shows that the blended [Fe II] emis- 
sion near 5200 A is weak, which indicates a relatively 
small 56 Ni mass (see $4] for more information on the 
mass of the ejecta). The line profiles of Mg I] A4571, 
[O I] AA6300, 6364, and [Ca II] AA7291, 7324 in the 
spectra of SN 2001ig in Figure [3] are quite strong and 
complex; Figure [4] illustrates them in detail. The Mg I] 
and [O I] lines at all epochs appear to have a broad 
peak, punctuated by multiple local maxima. Mg I] has 
peaks near —2000 km s _1 and at —1000 km s _1 rela- 
tive to the systemic ve locity of NGC 7424 (939 km s _1 ; 
iKoribalski et all [20041. This line may also have an- 
other bump near —3000 km s _1 in the blue wing of the 
main feature. The [O I] emission is even more complex, 
with local peaks near —3000 km s _1 , —1700 km s" 1 , 
150 km s _1 , and 1200 km s _1 (again relative to the sys- 
temic velocity of NGC 7424). The [Ca II] AA7291, 7323 
emission in Figure[4]has a broad peak which is blueshiftcd 
by ^1000 km s _1 in all three epochs. 

FigureQ]also shows the temporal evolution of the Mg I] 
A4571 and [O I] AA6300, 6364 line profiles. One of the 
most notable changes is the decrease in relative flux in 
the far wings of the Mg I] line near ±5000 km s _1 . As 
mentioned above, the feature near —3000 km s _1 appears 
to decrease in strength with time. This decrease of rel- 
ative flux with time near —3000 km s^" 1 is also seen in 
the [O I] feature. Finally, we note that there is also an 
increase in the flux in the wings of the [Ca II] feature 
near ±2000-5000 km s" 1 . 

The complex structure of the Mg I] A4571 and [O I] 
AA6300, 6364 lines has been seen previously in an- 
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Fig. 2. — Spectra of SN 2001ig (and a comparison SN lib) with some line identifications. From top to bottom, the spectra were obtained 
on 2002 Jan. 13, 17, 18 , and 29; ages shown are relative to the adopted explosion date of 2001 Dec. 3. The bottom spectrum is SN 1993J 
(Filippcnko ct al. 1994). Recession velocities have been removed from the spectra. Telluric features are indicated by the © symbol. 



other well-studied Type lib, SN 1993J (|Matheson et all 
120001) . They found that the [O I] and Mg I] line pro- 
files resulted from clumps in the ejecta while the [Ca II] 
lines did not. They also pointed out that this is con- 
sistent with t he explosion s cenar i o firs t proposed for 
SN 1987A by iLi fc McCravi lfl99l fl993h . in which the 
emission from oxygen comes from clumps of newly syn- 
thesized material while the calcium emission originates 
mainly in preexisting unifo rmly distributed material. 
Filippcnko & Sargent (1989) found similar irregularities 
in the [O I] AA6300, 6364 lines of the SN lb 1985F. They 
suggested that the spectral peculiarities come either from 
Raylcigh- Taylor fingers of higher-speed material or from 
local density enhancements. 

The [O I] doublet of SN 200 lig somewhat resembles 
that of the peculiar SN lb 2005bf (jModiaz et al.ll2008l) 
in that SN 200 lig may also show a double-peaked profile 
with one of the peaks near zero velocity and the other 
blueshifted by a few thousand km s _1 . However, the 
Mg I] line in SN 2001ig at all epochs has a broad profile 
with a complex top and does not seem to exhibit the same 
double-peaked natu re of the [O I] emi ssion. This differs 
from the results of iFolev et all (|2003t ) for the energetic 
broad-lined SN Ic 2002ap, which showed identical [O I] 
and Mg I] profiles. In addition, the flux in the Mg I] line 
appears to be more centrally concentrated in SN 2001ig 
compared to SN 2002ap. 

It has been proposed that the double-peaked [O I] pro- 
file of SN 2005bf w as caused by a unipolar blob or jet 
(|Maeda et al.H2007T ) while its polarization signature was 
attributed to a possible "tilted jet" passing through an 



asymmetric photosphere (|Maund et alj l200"7al lbl) . Both 
our models of late-time spectra of SN 200 lig and the 
similarity of th e profiles in our late-tim e spectra to those 
of SN 1993J (jMatheson et all I2000D indicate clumpy, 
oxygen-rich ejecta as well (see Sj4] for more information 
regarding the models and ejecta). Thus, SN 2001ig prob- 
ably has large-scale asymmetry in its ejected oxygen, ei- 
ther in the form of a blob or jet (tilted or otherwise), 
while the distributions of magnesium and calcium ap- 
pear more uniform (though the magnesium seems only 
marginally more uniform than th e oxyg en) . This is sup- 
ported by iFransson fe Chevalier! (|1989f ) , who predicted 
that magnesium and calcium should be found closer to 
t he core of the explosion t han ox ygen . 

iFilippenko fe Chornockl (|2002l ) stated that the Mg I] 
A4571 feature was the strongest (relative to [O I] AA6300, 
6364) ever seen in a SN. We investigate this claim more 
closely and the results are shown in Tabled where we list 
the relative line strengths for each nebular spectrum of 
SN 20 01ig along with some of the values from lFoley et alJ 
(|2003l Table 6) and values calculated from other SN 
spectra in the literature. To compute the relative line 
strengths for SN 2001ig we subtracted a local continuum 
from each spectral feature and then integrated the flux 
in each line . Thi s is the same procedure presented by 
IFolev et all ((20031) . 

Although it has been pointed out that this method 
of line-strength m easurement is somewhat inaccurate 
(|Folev et al.l 120031 ), it does seem to indicate that the 
Mg I] line of SN 2001ig grows with time with respect to 
the [Ca II] AA7291, 7324 doublet and the [O I] doublet. 
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Fig , 3. — Spectra of SN 2001ig (and two comparison SNe) with some line identifications. The top spectrum is SN 1993J (Filippcnk o et al.l 
119941) . The next three spectra are SN 2001ig and were obtained on 2002 Aug. 16 (from Maund et al. 2007b, Figure 2, reproduced with 
permission), Oct. 8, and Nov. 8; ages shown are relative to the adopted explosion date of 2001 Dec. 3. The bottom spectrum is the Type Ic 
SN 2002ap (Foley et al. 2003J). Recession velocities have been removed from the spectra (657 km s" 1 for SN 2002ap, from Lu ct al. 1993). 



TABLE 2 

Integrated Fluxes Relative to [Ca II] AA7291, 7324 



Supernova 


Day a Mg I] 


[O I] 


Ca II near-IR Mg/O Ratio 




4571 A 


6300, 6364 A 


8498, 8542, 8662 A 



SN 2001ig(IIb) b 


256 


0.38 


1.77 


C 


0.21 


SN 2001ig(IIb) 


309 


0.74 


2.06 


0.30 


0.36 


SN 2001ig(IIb) 


340 


0.83 


1.99 


0.57 


0.42 


SN 2008D (Ib) d 


109 


0.34 


1.02 


5.61 


0.33 


SN 2006aj (Ic-BL) C 


153 


1.19 


4.39 


1.37 


0.27 


SN 2003jd (Ic-BL) f 


367 


0.73 


1.58 


0.80 


0.46 


SN 1993J (IIb) s 


298 


0.27 


1.79 


0.36 


0.15 


SN 2002ap (Ic) h 


242 


0.72 


2.17 


0.91 


0.33 


SN 2002ap (Ic) h 


274 


0.63 


2.01 


0.66 


0.31 


SN 2002ap (Ic) h 


336 


1.00 


2.29 


0.45 


0.43 


SN 2002ap (Ic) h 


386 


0.87 


2.19 


0.46 


0.40 


SN 1998bw (Ic) h 


215 


0.23 


1.47 


0.41 


0.15 


SN 19941 (Ic) h 


146 


0.26 


0.97 


0.84 


0.26 


SN 1985F (Ib/c) h 


~240 


0.26 


3.15 


0.56 


0.08 



a Days since explosion.* 3 Spectrum previously published by [Maund et al. (2007b|). c Spec- 
trum ends at 8600 A. d Spectrum previously published by IModiaz et al.l ll2009T ), e Spec- 
trum previously published by Mazzali ct al. ( 2007a). ' Spectrum previously published by 
IValenti et a l. (2008). g Spectrum previously published by Filippcnk o et all lll994l) . h Values 
from lFolev et alj l|2003l. Table 6). 
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Velocity (1000 km s _1 ) 

Fig. 4.— The profiles of Mg I] A4571 (top), [O I] AA6300, 6364 (middle), and [Ca II] AA7291, 7324 (bottom) in SN 2001ig at three epochs, 
relative to the systemic velocity of the host galaxy NGC 7424. Ages relative to the adopted explosion date of 2001 Dec. 3 are shown. To 
create these figures, a local linear background was subtracted from each feature and the maxima of the features were all scaled to 1.0. 
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Relative to the [Ca II] feature, the Mg I] emission from 
SN 2001ig is somewhat weaker than that of SN 2002ap 
at similar epochs . The other SNe that also appeared in 
iFolev et all (|2003l . Table 6) all have about the same Mg I] 
strength, which is lower that that of either SNe 2001ig 
or 2002ap. However, we also perform this measurement 
on a few SNe from the l iterature (see Table d) and find 
that while SNe 2008D (jModiaz et al.ll2009ft and 1993J 
(|Matheson et al.ll2000ft have comparable Mg I] strength, 
this feature is quite strong in the broad-l ined Type Ic 

Ic-BL) SNe 2006a j ()Mazzali et all l2007al) and 2003jd 

Valenti et aT1l2008f ). 

Again comparing S N 2001ig to the SNe that appeared 
in lFolev et all (2003, Table 6), it has a Mg/O ratio com- 
parable to that of SNe 19941 and 2002ap, but larger than 
that of SNe 1998bw and 1985F. In addition, the Mg/O 
ratio of SN 200 lig is similar to that of SNe 2008D , 2006aj , 
and 2003jd (which in fact has quite a large ratio), but 
significantly greate r than that of SN 1993 J . This con- 
firms the claim of IFilippenko fc Chornockl (|2002f ) that 
SN 2001ig had the strongest nebular O/Mg emission-line 
ratio of any SN observed at that time. Table[2]also seems 
to indicate that it is not uncommon for core-collapse SNe 
to have relatively large Mg/O ratios, yet what distin- 
guishes between the SNe that do or do not show this is 
unclear. 

4. SPECTRAL MODELS 

By several months after the explosion of a SN, the den- 
sities in the ejecta are sufficiently low from expansion 
that the opacity drops below unity and the gas obeys 
nebular physics. At this time it is heated by radioactive 
decay of 56 Ni and cooled by line emission. The result- 
ing line profiles map rather directly the distribution of 
the heating material and can be used to derive several 
properties including 56 Ni mass, mass of t he emitting el- 
emen ts, and morphology of the explosion (Mazzali et al. 
120051) . 

In order to derive accurate models, we need our spectra 
to be as spectrophotometrically well calibrated as possi- 
ble. While the relative fluxes of our data are accurate 
because of our flux-calibration routines (see $2}, the ab- 
solute fluxes may be off due to imperfect observing con- 
ditions. To address this, we obtained publicly available 
photometric data on SN 2001ig from ESO and the Space 
Telescope - European Coordinating Facility (ST-ECF) 
Science Arc hive Facility 9 . The observations were single 
60 s Bessell (|Besselllll990T) i?-band images taken with the 
FORS1 spectrograph mounted on th e ESO 8.2-m Very- 
Large Telescope Melipal and Antu (|Ap pcnzcllc r et all 
[19981) . 

We calibrated both images, observed on 2002 Oct. 12 
and 2002 Dec. 8, to the United States Naval Observatory 
B1.0 (USNOB) Catalog with three comparison stars and 
measured magnitudes for SN 2001ig of 17.99±0.08 and 
19.13±0.14, respectively. In addition to our stated statis- 
tical uncertainty, we adopt a systematic uncertainty of ~ 
0.25 mag based on comparisons of USNOB magnitudes to 
Bessell i?-band magnitudes obtained a s part of the Lick 
Obse rvatory Supernova Search (LOSS; IFilippenk o et alJ 
1 2 If) photometry follow-up program. We then synthe- 
sized Bessell i?-band magnitudes from our late-time spec- 

9 http://archive.eso.org/ 



tra (taken on 2002 Oct. 8 and 2002 Nov. 8) and compared 
them to the magnitudes derived from the photometry, 
assuming a linear decline in magnitude during the last 
three months of 2002. In order for our spectra to match 
the photometry, and thus for the absolute flux to be accu- 
rate, our spectra from 2002 Oct. 8 and 2002 Nov. 8 must 
be scaled up by factors of 4.41 and 1.29, respectively. 
Due to the photometric measurements alone, the uncer- 
tainty in the absolute flux of each of these two spectra is 
-25%. 

We modeled both late-time spectra of SN 2001ig us- 
ing a code that computes nebular-line emission from a 
gaseous cloud. The gas is heated by the deposition of 
7-rays and positrons produced by the decay of radioac- 
tive 56 Ni (which is produced in the SN explosion itself) 
into 56 Co and then into stable 56 Fe. The energy thus 
deposited in the gas is thermalized collisionally, which 
leads to excitation and ionization of the gas. Heating 
is balanced by cooling, which takes place via line emis- 
sion. Most of the emission lines are forbidden, but some 
permitted transitions (e.g., the Ca II near-IR triplet) are 
als o active. The c ode is based on the original description 
bv lAxelrodl (|1980l ). The propagation o f the 7-rays is fol- 
lowed with a Monte Carlo procedure (jCappellaro et al.l 
1 19971 ) . Level populations are computed assuming non- 
local thermodynamic equilibriu m (non-LTE). A more d e- 
tailed description is provided bv iMazzali et al.1 (|2007bl ). 

Modeling the spectra also requires knowledge of the 
reddening of, and distance to, the SN. As discussed 
in SJSHI the reddening within NGC 7424 is negligible 
and we have only an upper limit from our data, so we 
adopt a total redde ning of E (B — V) — 0.011 mag from 
iSchlegel et al.1 (|1998l ) for the models. For the distance to 
SN 2001ig we ha ve used a recession velo city of 939 km s" 1 
(z = 0.003132; iKoribalski et~aT1 l2004h throughout this 
paper, yielding a distance mod ulus of m — M = 30.5 mag 
with H = 73 km s" 1 Mpc" 1 (jRiess et al.ll2005| ). assum- 
ing the entire recession velocity is due to the Hubble 
flow. However, there is some uncertainty in this value 
due to the relatively small distance to NGC 7424. A 
search of the literature reveals a range in distance moduli 
for this galaxy of 30.2-30.7 ma g using various indepen- 
dent t echniques (HyperLEDA 10 ■jTullvj 19881 : iMould et all 
120001 : iMeurer et all 120061 : lDavisll2009L private communi- 
cation 11 ). The value of 30.5 mag that we use here is a 
reasonable average as long as we keep in mind that it 
has an uncertainty of —0.3 mag which leads to an un- 
certainty of ~30% in the absolute flux of each late-time 
spectrum. When this is combined in quadrature with 
the uncertainty from the photometric measurements, the 
total uncertainty in the absolute flux of each of our late- 
time spectra is ~39%. 

The late-time spectra and models are shown in Fig- 
ure El both models appear to fairly accurately reproduce 
the major features observed in the data. However, the 
emission near 4000 A and the O I feature near 7700 A 
seem somewhat suppressed in the models. In addition, 
the intensity of the Ca II feature near 8600 A in the mod- 
els is somewhat larger than that of the observed spec- 
trum, leading to a small inconsistency in the position of 

10 http://leda.univ-lyonl.fr/ 

11 Using the 2MASS Redshift Survey (2MRS) galaxy catalog to 
model bulk flows in the local Universe. 
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the combined Ca II plus [C I] A8727 emission. 

Here we use the one-zone ver sion of the model de- 
scribed by [Maiiair]i~ai] ((20073) • It requires an emit- 
ting nebula with an outer velocity of 4300 km s _1 , as 
determined from reproducing the width of the strongest 
emission lines ([O I], Mg I], and [Ca II]). The spectrum is 
powered by —0.13 Mq of 56 Ni. This is determined from 
the combination of the intensity of the [Fe II] emission 
feature near 5200 A and the requirement that the energy 
emitted in the radioactive decay chain of 56 Ni excites 
other ions and reproduces the observed emission in other 
lines. This is a slightl y larger value than what was de- 
rived for SN 1993J (see lSwartz et al.lll993l . and references 
therein), but it is still fairly typical for core-collapse SNe. 

The total mass enclosed within the outer velocity of 
4300 km s _1 is ~ 1.15 M . Interestingly, our models of 
SN 2001ig are quite similar to that of the SN Ic 19941 
which underwent significant mass loss, most likely in a 
binary interaction, and w h ich also ejected ~lM a of ma- 
terial (jSauer et al.ll2006h . iNomoto et "all (|1994l ) showed 
the progenitor star of SN 19941 to be — 15 Mq initially, 
which suggests that the progenitor of SN 2001ig was also 
initially -15 Mr. 

Furthermore. I Van Dvk et all (|2002t ) and lMaund et all 
(|2004h found that both the companion and the progeni- 
tor of SN 1993J had intial masses of —15 Mq and that the 
companion of SN 1993 J is currently 22 M©. iRvder et al.1 
(12006ft discovered that the companion of SN 200 lig is 
currently 10-18 M Q , so it seems likely that the pro- 
genitors of SN 200 lig and SN 1993 J (and possibly even 
SN 19941) had similar initial masses, and that their com- 
panions have similar masses as well. The differences in 
the companions' masses (as well as physical separations 
from the progenitors) probably led to slightly different 
mass-loss and interaction histories which gave rise to the 
spectroscopic differences observed. 

Table [3] shows the complete mass composition of the 
models. The upper and lower error bars represent the 
change in the derived abundance after increasing and 
decreasing the measured flux by 39% to reflect the flux- 
calibration uncertainties mentioned above. As in the 
nebular spectra of all other stripped-envelope SNe, and 
unlike the case in SNe la, [Fe III] lines are missing (in 
particular a strong feature near 4800 A). This is best re- 
produced by assuming fairly significant clumping in the 
eject a, which is a signatur e of an aspherical expl osion 
(e.g.. iMazzali et all [20070 ). iMaund et all (|2007bl ) use 
spectropolarimetric observations of SN 200 lig to come 
to the same conclusion regarding the asphericity. 

The dominant element in the ejecta by mass is oxygen 
(0.81 Mq). The derived carbon abundance is quite small 
(0.04 Mq); otherwise, the emission near 8500 A would 
become too strong. We point out that a small C/O ratio 
is not unusual in stripped-envelope SNe. In addition, 
hydrogen may be present only marginally at these low 
velocities, as deduced from the almost complete absence 
of any emission (as seen in Figure [3|). 

As mentioned in §3.31 the late-time spectra of 
SN 2001ig have an unusually large Mg/O intensity ratio, 
comparable to that of SN 2002ap. However, the Mg/O 
ratio by mass of SN 2001ig is -0.014 (see TableEB while 
that of SN 2002ap is -0.032 (|Mazzali et al.ll2007bl . This 
difference in masses is consistent with the fact that the 



TABLE 3 
Mass Composition of the Models 51 



Element 


Mass (M Q ) 


Mass (M ) 




Day 309 


Day 340 



c 
o 

Na 
Mg 
Si 
S 

Ca 

56 Ni 

Total 



2.0t°;g X 10" 2 
8.1^2.? x 10 _1 
7.0±5'5 x !0~ 5 

x to-* 

1.0t°-° x lO" 1 
3.0±0-g x lO" 2 

3.1+o 9 x !0~ 2 
1-5+gjxlO-i 

1 n + U.a 

L - z -0.3 



5.0+.°;° x io- 2 

8.2t\^ x IO" 1 
1.0±g'g x 10" 4 

i-i±8:S >< i°- 2 

l.OlH x IO" 1 
3.0+1° x 10-2 
4.Qto 9 x 10~ 2 
l-l^xlO- 1 

1 9 +U.a 

L -^-0.2 



a The upper and lower error bars represent the 
change in the derived abundance after increas- 
ing and decreasing the measured flux by 39% to 
reflect the flux-calibration uncertainties men- 
tioned in i|4l 

envelope of SN 2002ap was extremely stripped before ex- 
plosion, since it s howed no spectros copic evidence of hy- 
drogen or helium (|Folev et al.ll2003T ). On the other hand, 
the envelope of SN 2001ig was only partially stripped 
prior to core collapse since we observed a significant 
amount of hydrogen and helium at early times. 

5. CONCLUSIONS 

In this paper we have presented and analyzed optical 
spectra of SN 2001ig. One week after explosion, the SN 
defied simple spectral classification, but by two weeks 
after explosion SN 200 lig appeared to be part of the 
Type lib subclass of SNe. This was confirmed as the 
SN went through a transition from Type II to Type lb 
during the first few months of observations. 

By nine months after explosion, SN 2001ig had entered 
the nebular phase, revealing some of the strongest Mg I] 
A4571 and [O I] AA6300, 6364 features ever observed in 
a SN; its Mg/O intensity ratio is one of the largest ever 
seen as well. 

We derive models of SN 2001ig from our spectra taken 
309 and 340 days after explosion, showing that the ma- 
jority of the inner ejecta (below 4300 km s^ 1 ) were in 
the form of oxygen and a significant fraction (nearly 
0.13 M Q ) in the form of 56 Ni. Additionally, there ap- 
pears to be a distinct lack of hydrogen in the inner ejecta. 
We also find that the total mass of this inner ejecta was 
— 1.15 M , suggesting that the progenitor of SN 2001ig 
was a relatively low-mass star (—15 Mq). 
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Fig. 5.— The model and observed spectra of SN 2001ig from 2002 Oct. 8 (309 days after explosion, top) and 2002 Nov. 8 (340 days 
after explosion, bottom). The recession velocity of the host galaxy NGC 7427 has been removed from the spectra. The model accurately 
reproduces the major spectral features seen in the data. 
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